2x 
“JTA —1 
Walaupun bentuk analitik fungsi pemeta z(C) tidak diketahui tetapi dapat ditaksir oleh 


OL, -(j-1) 





(4.124) 


sebuah fungsi kubik lentur dalam variabel kompleks, yaitu fungsi W(C). Persyaratan 


kesinambungan turunan kedua dari WC), yaitu nilai W'' diberikan oleh persyaratan 
(4.109) dalam variabel kompleks, yaitu untuk 2 £ j  JTA —1 
jl jl 


itna an asaa,.-sa, nat 
| j j 





AGL, WA ACL,  # ACL, ). WI 4 ACL,. WI, 6 aa z 2 nara 
(4.125) 
Nilai W, adalah nilai ZAA , yang diberikan, yaitu 
W, -ZAA, (4.126) 
sedangkan AGL, diberikan oleh rumus berikut 
AG 5 GL —GLy (4.127) 


mengingat bahwa untuk aerofoil berekor lancip titik ekor adalah titik kritis, yaitu 
2/6) — 0, maka 

W - Wi, 50 (4.128) 
Ini berarti bahwa kondisi batas yang berlaku adalah untuk kasus II dan sistem persamaan 
(4.125) kemudian dapat diselesakan seperti telah dijelaskan sebelumnya, tetapi dengan 
mengingat bahwa koefisien aj, bj, c, dan d, adalah variabel kompleks. 


Rumus interpolasi untuk selang j diberikan oleh rumus (4.121) dalam variabel 
kompleks, yaitu 


war) - (e1,.(au —tu,) 02) au, A03). —CL,)-ZAA, 
(4.129) 
dimana koefisien kompleks C1, C2 dan C3 adalah 
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co, Waw 
“6 ACL, 
1 
C2 Wi 
C3, - W, 


Panjang kurva sepanjang permukaan aerofoil dari titik ZAA, ke ZAA ,,, diberikan oleh 


rumus eksak berikut 


jl jl 


AS, -SAA,,, -SAA, — | ds— ag 
i j 
dz 


AS, — ja ace 


yang ditaksir oleh rumus eksak berikut 


ac, 


aer 0 - Jiekao 














AS, - iw'tenjaor (4.130) 
j 


selang AOL — 21/(JTA —1) antara OL, dan OL,,, kemudian dibagi menjadi 2N sub- 
selang yang lebih kecil sebesar 50 L, dimana 80 L — A0L/2N. 

Integral (4.130) kemudian dapat dihitung dengan menggunakan rumus Simpson 1/3 
sebagai berikut 


Ra 
AS, - an (s “ E, #2), Ea “su 


dimana 


x -Iwi|-Iw(an,)|-Iw'expliox,)) 





Iwa) -Iwepliou,,.)) 








Ea 5 (Wa 
| Anta Iw'(explilox, k01) 


won) - (sle1,-(au—an,) 202). 03.) 
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Jarak sepanjang permukaan aerofoil dari titik j-1 ke titik j kemudian dapat dihitung 
dengan rumus 

SAA, —0 

SAA,,, -SAA, AS, untuk 1S j« JTA-1 
Setelah distribusi panjang permukaan aerofoil SAA, diperoleh kemudian kita dapat 


menggunakan himpunan pasangan nilai (xaa,,saA,), untuk 1S jSX JTA, untuk 
menghitung koefisien-koefisien fungsi interpolasi kubik lentur x(s). Disamping itu 
koefisien-koefisien fungsi interpolasi kubik lentur y(s) yang melewati titik-titik data 
(vaa,,sAA,) juga dapat dihitung. Fungsi-fungsi interpolasi x(s) dan y(s) kemudian 
dapat digunakan untuk menghitung koordinat titik-titik sembarang pada permukaan 
aerofoil, yaitu (XAI,YAI), adalah jarak sepanjang permukaan aerofoil dari titik ekor 
ZAA ke titik interpolasi, ZAI, diberikan sebagai berikut 


XAI — x(s — SAI) 
YAI -— y(s — SAI) 


Fungsi interpolasi x(s) dan y(s) adalah polinom kubik lentur dalam variabel riil dan cara 
menghitung koefisien-koefisiennya telah dijelaskan sebelumnya. 
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BAB V 
KOEFISIEN PENGARUH PANEL SINGGULARITAS 


Dalam bab-bab sebelumnya kita telah membahas masalah aliran potensial, tak 
mampat dengan menggunakan bantuan teori variabel kompleks dan pemetaan konformal. 
Persamaan atur atau model matematis masalah aliran potensial tak mampat adalah 
persamaan Laplace yang bersifat linier. Sifat linieritas ini dapat dimanfaatkan untuk 
memperoleh solusi yang rumit sebagai kombinasi linier dari solusi-solusi sederhana atau 
elementer yang diketahui rumus analitisnya. Solusi persamaan Laplace adalah semua 
fungsi analitik komplek yang dapat diturunkan dengan bantuan teori variabel kompleks. 
Beberapa contoh dari fungsi analitik kompleks elementer yang bermanfaat dalam 
aerodinamika adalah aliran sumber/serap, doblet dan vorteks yang telah dibahas secara 
rinci di modul I buku bahan kursus ini. Fungsi potensial kompleks untuk sumber, vorteks 
dan doblet bersifat analitik di semua titik kecuali titik dimana solusi elementer tersebut 
berada. Titik dimana fungsi potensial kompleks menjadi tak-terdefinisi disebut titik 
singguler atau singgularitas dari fungsi tersebut. Oleh karena itu sumber, doblet dan 
vorteks, yang telah dibahas adalah jenis diskrit atau cacah karena diwakili oleh titik-titik 
diskrit, dan oleh karena itu disebut singgularitas diskrit. Dalam bab ini kita akan 
membahas tentang singgularitas terdistribusi kontinyu (mulus) atau singkatnya 
singgularitas kontinyu yang merupakan konsep lanjutan dari singgularitas diskrit. 
Seandainya titik P dan O kita " letakan " singgularitas sumber, masing-masing dengan 
kekuatan tP dan 1O, maka fungsi potensial kompleks di titik T dengan koordinat ZT 
adalah : 


k10) 
In 

2m 

dimana ZP dan Z0 adalah koordinat dari titik P dan 9. Pada umumnya apabila ada 





01) — In(zT ZP)1 Sin(zT-Z0) (5.1) 


sejumlah N sumber, masing-masing berada di titik S(n) dengan kekuatan 1S(n), untuk 
ISn£N, maka fungsi potensial kompleks di T adalah 


&(T)- - 1S(n).in(zT —ZzS(n)) (5.2) 


dimana ZS(n) adalah koordinat titik S(n). Sekarang marilah kita bayangkan kasus dimana 
semua titik S(n) berada pada garis yang menghubungkan titik P dan titik 9. Garis PO 
tersebut dibagi menjadi N selang seragam dengan panjang selang sebesar 


As -POy (53) 
Titik P adalah salah satu ujung selang pertama dan 9 adalah ujung selang ke-N atau 


selang terakhir. Titik tengah selang ke-n adalah titik S(n) sedangkan s adalah jarak dari 
titik P ke titik S sepanjang garis PO. Jarak dari titik S(n) ke titik T diberi simbol r(n) 
sedangkan @(n) adalah sudut yang dibentuk oleh garis S(n)T dengan garis PO yang 
sejajar sumbu-x, jadi 

ZT—ZS(n) —r(n).exp(i0(n)) (5.4) 


Rumus (5.2) untuk fungsi potensial kompleks sekarang dapat ditulis ulang menjadi 
N TS 
D(T,PO) - - Sa r(n)--i0(n)).As (5.5) 
apabila nilai N diperbesar, diambil limitnya mendekati tak berhingga maka As akan 


n-l 


1S(n) 





mendekati nol, tetapi nilai ditentukan bernilai tetap berhingga, yaitu 





PO limit Sa Ap) (5.6) 


As—0 S 
Perumusan (5.6) menunjukkan bahwa t(n) adalah kekuatan sumber per satuan panjang 


di titik S(n) dan mempunyai nilai tertentu untuk masing-masing n. dengan pengertian dan 
perjanjian tersebut persamaan (5.5) sekarang dapat ditulis menjadi 


—PO 
Kn at , 
OPO) limit Da) lite) -i0ln)).As (5.7) 
Dari teori kalkulus dasar dapat diketahui bahwa persamaan (5.7) dapat diganti oleh 
integral berikut 
(1,20) - | o)im3) sio) 68) 
T 


P 
Penalaran yang sama dapat juga diterapkan untuk singgularitas vorteks dengan hasil 
berikut 


kel 


&(T,P0) --— | 10) ilnr(s) 0(s)).as (5.9) 


P 
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Sedangkan untuk doblet rumus yang diperoleh adalah 


O(T,PO) z 5) K)-oy (00506) —isin @(s))-as (5.10) 


dimana y(s) dan y(s) masing-masing adalah kekuatan vorteks dan doblet per satuan 


panjang di titik S. Dalam rumus di atas sudut @ bernilai positif untuk titik T yang berada 
di atas garis PO dan negatif bila T berada di bawah garis PO, dengan besaran O sampai 


dengan 1, atau 





o s/o(s)| xx (5.11) 
yaitu 

OxO(s)«T untuk T di atas PO 

026(s) - untuk T di bawah PO 


Komponen-komponen real dan imajiner dari @(z) adalah fungsi potensial &(r,@) dan 
fungsi arus w(r,@) yaitu 


O(T,PO) - H(T,PO) -w(T,PO) (5.12) 
Untuk panel P9 yang cukup pendek, distribusi kekutatan singgularitas per satuan 
panjang k(s) dapat ditaksir secara akurat oleh polinom berderajat rendah sebagai berikut 


kls) - kk 50) t13( K0) HOT (5.13) 


dimana HOT adalah "higher Order Terms" atau suku-suku deret berderajat lebih tinggi. 
Dengan pengertian tersebut rumus-rumus untuk fungsi potensial dan fungsi arus 
yang diimbas oleh panel singgularitas kontinyu dapat ditulis sebagai berikut 


Sumber/serap: 


6(T,PO) — Zlsho r.ds ea Kajar.as 01) 


w(T,PO) -— Lose | oko Jo euor| 
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Vorteks: 


6(T,PO) - £- -njo. al Yo pas euor!| 


w(T,PO) — &- snfmeas- k0)araerwor| 


Doblet: 





H(T,PO) - Fla Pawan 00) 
w(T,PO) 2 aja s0) 


Rumus-rumus di atas dapat ditulis ringkas sebagai berikut 








Sumber/serap: 


6(T,PO) — Yth) 4#HOT 


k-i 


wT,PO)-— Ha AN O,k— 1) -HOT 


k-1 
vorteks: 


K 
6(T,PO) ja —yk.X2,k —1) #HOT 


21 


K 
w(T,PO) - jk. Kt,k—1) HOT 
T 


k-i 
Doblet: 


6(T,PO) -— sa Aa aan 1) -HOT 


k-1 


K 
w(T,PO) jay —pk.1(4,k — 1) #HOT 


21 
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Jana 


(5.14) 


(5.15) 


(5.16) 


(5.17) 


(5.18) 


(5.19) 


dimana 


X1,k—1) Noka) Inr.ds (5.20) 
X2,k—1) Nota) ca (5.21) 


k-1 


V3,k — »- X0) (c080/ ).as (5.22) 
W4,k —1) — Iko y (sing/).as (5.23) 


Disamping informasi tentang fungsi potensial dan fungsi arus, dalam masalah 
aerodinamika biasanya dibutuhkan juga informasi tentang distribusi kecepatan di medan 
aliran. Vektor kecepatan itu sama dengan konjugat dari turunan fungsi potensial 
kompleks, yaitu 
3 aan 
VT)- ad - @'(T) 
dz 
V(T) -u(T) -i.v(T) 
dimana u dan v masing-masing adalah komponen kecepatan searah dengan PO dan 
tegak lurus (normal) pada PO. Konjugat kompleks dari turunan fungsi potensial 


kompleks adalah 
Sumber/serap: 














rara d ( T nz) T 1 exp(-i0) .: exp(i0) 
dz 27 211 21x r 21 r (5.24) 
3 1 cos@ .1 sing 
U-tiv ——1. Ai —T 





2x € “2x” r 
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Vorteks: 


aa) ane) ng —iy exp(-i0) y .xp(i0) 











dz 2x 21 r 21 r (5.25) 
2 1 sin@d .1 cos 
u-iv— : Hi. s 
2m r 2x r 
Doblet: 





“dal 2m Tz” 2x r? 27 (5.26) 


1 cos290 1 sin20 


2 


&'z) an u exp(-i20) jp exp(i20) 





u-—iv 





PN Dada” 
Untuk singgularitas kontinyu rumus-rumus di atas dapat diolah menjadi 


9 
(056 Jas 4 | Il :000 as t or! 
r SPO r 


V(T,PO) Ha anjar 25 ons) nor| 


Vorteks: 


u(T,PO) 5 ni Anas in 


P P 


Oo 
V(T,PO) - Hn ejaan 2 e2)asenor| 


Doblet 


ko) kol 
u(r,po) 2 (oran pa 2)acrnor!| 
P P 


Sumber/serap: 





u(T,PO) tak 


“——0 
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V(T,PO) - : ai ino), ds si ja ds #01! 


P 


Rumus-rumus di atas dapat di tulis ringkas sebagai berikut 








Sumber/serap: 
u(T,PO) -— ai perak 1) -HOT (5.27) 
V(T,PO) — 5 ae M2 1) -HOT (5.28) 
T x-1 
Vorteks: 
K 
u(T,PO) - LH Y(4,k —1) -HOT (5.29) 
T x-1 
V(T,PO) -— Ta yk. 3, — 1) -HOT (5.30) 
k-1 
Doblet: 
1 K 
u(T,PO) -— Y —uk.YS,k—1) #HOT (5.31) 
21 
1 K 
VT,PO) —-— X tk. I(6,k —1) # HOT (5.32) 
k-1 
dimana 
“(ss sp cos20 
XU5,k—1)- (| —| . .d 5.33 
Ga) Ce ia 
8( s ea sin 20 
W6,k—1)- (| —| . .d 5.34 
G0) Aja 2 


Integral I(l,k-1), (2,k-1) dan seterusnya dikenal sebagai koefisien pengaruh panel P9 
pada titik T dan dapat diolah lebih lanjut sebagai berikut. 
Perlu diingat bahwa s adalah jarak dari P ke titik S pada garis PO dan bernilai 


positif arah PO. Variabel r adalah jarak dari S ke titik T dan selalu bernilai positif. Sudut 
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» 4 
0 adalah sudut yang dibentuk oleh garis ST dengan arah PO dan nilainya memenuhi 
persayaratan (5.11). Untuk lebih memperjelas masalah, hubungan geometris antara titik- 


titik dan variabel-variabel yang terlibat ditunjukkan pada gambar (5.1). Garis PO 
ditunjukkan sejajar arah sumbu-x positif dan T berada di atas PO, yaitu vektor 
5» 5» , 
PO x PT mengarah keluar dari halaman buku ini menuju ke arah mata pembaca, 
sedangkan besarnya adalah 
PO x PT-PO.PT.sinTPO 
AN A 

dimana TPO adalah sudut positif, 0 c TPO « x. 

Kondisi yang diberikan ditunjukkan pada gambar (5.1a). Pernyataan "di atas" dan 


"di bawah" itu relatif, karena untuk kasus PO searah dengan sumbu-x negatif (lihat 
gambar (5.1) titik T yang kelihatannya berada di bawah garis P9 secara matematis 
berada di atas garis PO, karena untuk sudut 0 yang memenuhi persyaratan (5.11) vektor 


PO x PT bernilai positif. Sebaliknya pada gambar (5.1d) titik T kelihatannya berada di 
atas PO, tetapi sesungguhnya berada di bawah PO, karena untuk kasus ini vektor 


—» -D 
PO x PT bernilai negatif. Jadi sesungguhnya perkataan T berada di atas atau di bawah 
PO itu ditentukan oleh kenyataan apakah TR bernilai positif atau negatif, dimana R 


adalah proyeksi T pada garis PO sehingga 


dinr — PRE, 2 
“NG « 





dimana simbol caping “ berarti satuan vektor. 
Apabila koordinat titik-titik P, O dan T adalah (XP, YP), (XO,YO) dan (XT,YT) 


maka 


TR — ((XO-XP).(YT- YP)-(YO- YP).(XT-—xP))/PO 





atau diringkas menjadi 
TR — (AXOP.AYTP — AYOP.AXTP)/PO (5.35) 


dimana 
AXOP —- XO-— XP 


dan seterusnya. 

-» 5 » 
Panjang garis PR bernilai positif bila PR searah dengan PO dan bernilai negatif bila PR 
kebalikan arah PO. Nilai PR dengan demikian diberikan oleh rumus berikut 
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- PO.PT 
PR -?PO 37 


atau 


ja (AXOP.AXTP 4 aa Ka 
5 bg 


(5.36) 
Sekarang perhatikan bahwa titik S dengan koordinat (XS,YS) boleh berada di sepanjang 


PO dan s adalah jarak PS searah PO. 
Ini berarti bahwa 


XS — XP 4 (AXOP/PO).s 
YS — YP 4 (AYOP/PO).s 
Variabel r adalah jarak dari S ke T, yaitu 
r —(XT-xS) 4(YT- yS) 
r? - (XT- XP — (AXOP/PO).s) 4 (YT-— YP —(AYOP/PO).s) 


r? — (AXTP — (AXOP/PO).s) # (AYTP — (AYOP/PO).s)” 


r'—- PT'-—2PR.s 45? (5.37) 
dimana 

PT? — AXTP? # AYTP? (5.38) 
dari gambar (5.1) diketahui pula bahwa 

PT? — PR? 4 TR? (5.39) 
jadi rumus (5.37) dapat diolah menjadi 

r'- TR? 4(PR -s) (5.40) 
Panjang panel PO dan jarak OT diberikan oleh rumus-rumus berikut 

PO? — AXOP? # AYOP? (5.41) 

OT? — AXTO' 4 AYTO (5.42) 


Variabel OR diberikan oleh rumus 
OR -PO.OT/PO 
OR — (AXOP.AXTO # AYOP.AYTO)/PO (5.43) 


dan juga oleh rumus berikut 
OR-PO-PO (5.44) 


Selanjutnya perhatikan juga bahwa 
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sin — TR/r (5.45) 


cos0 — (PR —s)/r (5.46) 
tan0 — TR/(PR —s) (5.47) 
cos20 —1-— 2sin'@-—1—2TR”/r? (5.48) 
cos20 — 2sin0cos@- 2TR(PR —s)/r” (5.49) 


Rumus (5.20), (5.21), (5.22), (5.23), (5.33) dan (5.34) untuk koefisien pengaruh dapat 
diolah lebih lanjut sebagai berikut 




















p0-tattx-—0 4 Falestar)- Jesamr| (5.50) 
Pp. P 
PO .I2 &--t) dls or)- | | (5.51) 
12, -5 : $ 
P P 
9 k-1 k 
PO .1(3,k—1) — | PR . 4 — (5.52) 
Pp Pp 
&i s ds 
PO. 4,k—1)— | TR 5 (5.53) 
Pp 
Ta 8? ss , ds 
PO.I(S,k—1) J aa Ya Sa (5.54) 
P P 
0g k-1 k 
PO L.I(6,k —1) — | 2PR.TRS 25 Pangsa (5.55) 
Pp P 


Turunan dari persamaan (5.47) adalah 


dtan0 — d IR ) 





PR —s 
TR.ds 
1-#- tan? 9)/d0 - ——7 
( ) (PR —s) 





(ae) “Am 
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(R2 «(pr —5)')ao — TR.as 


ds 
OP STR (5.56) 


Rumus (5.40) dapat dimanfaatkan untuk menurunkan persamaan berikut 


1dr? 1d(PR—s) 



































di 
“2 r 
din PR SS (5.57) 
Selanjutnya perhatikan bahwa 
ds —dPR-s) 1 dcotd 
r' (TR24(PR—5)) TR (iteoro) 
ds 1 dg sin”0d9 
r“ TR'14cot@  TR' 
1—cos20)d0 
al 2 2) -— |ao-Lasin2e 
r 2TR 2TR 2 
ds d9 1 PR —s 
5.58 
rs 2TR' 2TR' ( r? ) (5-58) 
Persamaan (5.56) dan (5.57) dapat ditulis ulang menjadi 
ds d9 
TE — 5.59 
2 TR (5.59) 
aa) (5.60) 
r TR 


Kedua rumus di atas dapat digunakan untuk menurunkan rumus s"ds/r? untuk k 2 2 


secara beruntun sebagai berikut 


gas (P-PT'42PR.s).ds 





r r? 


2 
ka ds— PT? 05 opa Ss 
r r r 
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2 2 mp2 
S Os st 2PR-dinr#PR TR EN: 
r 





d9 (5.61) 


Untuk k-3 hasil yang diperoleh adalah 





3 2 

2 2 - s.ds—PT? 2 42PRS ba 
r r r 
s'ds 


SA — ds #2PR.ds # (3PR? TR Jan 4-3 (PR? — 3TR”).40 (5.62) 


2 

r 
Persamaan (5.57) dapat digunakan mengevaluasi suku ke-2 ruas kanan (5.50) sebagai 
berikut 





2 
sdinr ——PRSOS 4 59 
r r 
s.dinr — dst PR.dinr— TR.d9 (5.63) 
Untuk k-2 hasilnya adalah 
#dlnr — La #PR.ds-#(PR?—TR?)dinr— 2PR.TR.d0 (5.64) 
Suku ke-2 ruas kanan persamaan (5.51) dapat dievaluasi dengan bantuan rumus (5.56) 
sebagai berikut 
sds 
sd9 - TR—- - TR.dinr #PR.d9 (5.65) 
r 


Untuk k-2 hasilnya adalah 


s'ds 


r? 





Sd0 - TR“ - TR.ds #2PR.TR.dinr (PR? —TR?).40 (5.66) 


Persamaan(5.47) dapat digunakan untuk merumuskan sudut-sudut berikut 


00 —tan"(TR/OR) (5.67) 

OP —tan"(TR/PR) (5.68) 

A-00-0P (5.69) 
Mengingat bahwa 


tan(0O) — tan(OP) 
1-#-tan(0O).tan(OP) 
Maka rumus singkat untuk A dapat diperoleh sebagai berikut 





tan(0O — OP) — 
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TR(PR-OR) — PO.TR 









































tanA 5 5 
PR.OR4#-TR- PT —PO.PR 
yaitu 
A-tan” PIR (5.70) 
PT —PO.PR 
Rumus (5.58) dapat dikalikan dengan s dan hasilnya adalah sebagai berikut 
sds ”” —s.dlPR—s) — (-PR4(PR-s))a(PR -s) 
- "Ma 2 
Po (TR24(PR—s)) (rr2 «(pr —s)) 
d(PR — d(PR —s)” 
ng PR. ( ai ( s) - 
r (rr? -(2R—s)) 2 (rr? (pR—s)) 
sds PR-ds | Ldr” 
rs rs 2 rs 
s.ds PR 1 (1 PR PR —s 
— .dO——d d 557 
r—— 2TR? 2 (3) 2TR? | r? ) AN) 
Perkalian persamaan (5.58) dengan s? dapat memberikan hasil berikut 
2 r'—PT'42PR.s).ds 
s das | : ) La PT? 9S KPR dn 
r r r r r 
2 
s-ds d0 (2pR?—pr) 3s PR. : 
r TR r r 
hp 5 2 TR? ( PR— 
Sana pn TP Ge end (5.72) 
r 2TR 2TR r r 


Integral di ruas kanan persamaan (5.50) sampai dengan (5.55) sekarang dapat dihitung 
sebagai berikut : 
Untuk k-1 kita memerlukan rumus integral berikut 


0 
Ja(sinr) — PO.In OT 
P 
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te) 
Js-dinr - PO -PR.In(OT/PT)-TR.A 


Jatso) - P9.O9 


P 


kol 
|s-d0 — TR.In(OT/PT) #PR.A 
. 











jan 
be TR 

ko) 

(38 tm af Ju ERA 

br PT) TR 
jas 1 (OR PR 
br” 2TR' 2TR'(OT? PT 





pet PO.PR 
dr 2TR?' 2TR?.0T? 


Untuk k-2 kita membutuhkan beberapa rumus tambahan sebagai berikut 


0 
Jals.tnr) —pO?.mOT 


0 
| sdinr — OPO 42PR) 4 (PR? — TR?)in(OT/PT) -2PR.TR.A 
id 


9 
Jals:o) -po?.o9 


P 


9 
| tao - PO.TR 4 PR.TR.In(OT/PT) #(pr2 . TR?).A 


P 


9 2 Dm 2 
J IS P9 42PR.In(OT/PT) PR TR LE AN 
Ir TR 
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Je PT.A PO PO.OR 
! rt 2TR' 2TR?' 2TR?.OT 
Pemasukan hasil-hasil di atas ke-rumus-rumus (5.50) sampai dengan rumus (5.55) 
akhirnya memberikan hasil berikut : 
Untuk koefisien pengaruh panel PO di titik T. 
































Untuk k-1 
I(1,0) —-PO -4#PR.InPT-— OR.InOT-4# TR. A (5.73) 
I(2,0) — —TR In(OT/PT) -PR.OP -— OR.00 (5.74) 
(3,0 Na (5.75) 
X4,0) — (5.76) 
OR PR 
US,0) — Ara PT (5.77) 
TR TR 
16,0) — SP (5.78) 
Untuk k-2 
14,1) PO 2PR PO PRA TR Ore PR TR pp PRTR 
4 2PO PO 
(5.79) 
2 2 2 2 2 
12,1) TR Pa Ia) P3 — PR? TR 001PR -IR' pp 
2 PO PT 2PO 2PO 
(5.80) 
W3,1) ——1 5. af SE) EA (5.81) 
PO (PT) PO 
OT) PR 
(4,1) — ES 5 ia (gg 10 .A (5.82) 
X5,1) — 2 seb OT) (5.83) 
OT PO (PT 


3 


(5.84) 





16,1) - Aa Pan 


PO POL PT? OT 
Rumus-rumus di atas berlaku umum untuk semua titik T kecuali bila T berada 
pada perpanjangan garis PO, atau berada di antara titik T dan 9 sepanjang garis PO. 
Apabila T berada pada garis PO maka R berhimpit dengan T dan TR -0. Namun 
demikian PR tidak sama dengan PT karena PT selalu bernilai positif sedangkan PR 
boleh bernilai negatif walaupun nilai mutlaknya sama dengan PT. Oleh karena itu 
PR tetap diberikan oleh rumus (5.36) sedangkan PT diberikan oleh nilai mutlak PR yaitu 


PT — (PR) (5-85) 
Penalaran yang sama juga berlaku untuk OT yaitu 
OT - (OR) (5.86) 


dimana OR diberikan oleh rumus (5.43) atau (5.44). Sudut @0 dan OP diberikan oleh 
rumus (5.67) dan (5.68) dimana TR bernilai nol dan oleh karena itu 


00-—0 bila OR bernilai positif (587 
00-17 bila OR bernilai negatif - 
0P-—0 bila PR bernilai positif (5.88) 
OP-7x bila PR bernilai negatif i 


Nilai sudut A tetap diberikan oleh rumus (5.69), jadi 


A—0 bila PR dan OR bertanda sama 
An bila OR negatif dan PR positif (5.89) 
A-—-1 bila OR positif dan PR negatif 
Dengan pengertian tersebut sekarang persamaan (5.73) sampai dengan (5.78) dapat 
disederhanakan untuk kasus titik T berada pada garis PO menjadi 
Y1,0) ——PO 4 PR In|PR|— OR.In|OR| (5.90) 


5» 5» 
I(2,0)-0 bila OT searah dengan PO atau T disisi "kanan" O 
-1.OT bila T berada diantara P dan 0 (5.91) 


4 4 
-1.PO bila PT kebalikan arah PO atau T disisi "kiri" P 
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X3,0) ——In(|OR//IPR)) (5.92) 





I(4,0) —0 (5.93) 
KA PO 

X(5,0) — An AOR (5.94) 

I(6,0) —0 (5.95) 


Untuk kasus dimana titik T berhimpit dengan titik O, rumus-rumus di atas 


disederhanakan menjadi 


I1,0) —-PO--PO.InPO — PO(InPO 1) (5.96) 
(2,0) — 1.PO (5.97) 
I(3,0) — tak- terdefinisi (5.98) 
(4,0) — 0 (5.99) 
(5,0) -— tak- terdefinisi (5.100) 
I(6,0)-— 0 (5.101) 
Apabila T berhimpit dengan P rumus-rumus di atas menjadi 
T4,0) —-PO #PO.InPO - PO(InPO —1) (5.102) 
2,0) - 1.PO (5.103) 
(3,0) — tak - terdefinisi (5.104) 
X4,0)— 0 (5.105) 
W(5,0)-— tak - terdefinisi (5.106) 
I(6,0)-— 0 (5.107) 


Rumus-rumus untuk koefisien pengaruh kedua, yaitu persamaan (5.79) sampai dengan 
(5.84) juga dapat disederhanakan untuk kasus T berada disepanjang garis PO menjadi 
POt2PR  PO'-PR 


Y4,1)- - At EaO JnfOR| 175  tnfPRI (5.108) 








12, 1) -0 bila OT searah dengan PO atau T di sisi kanan O 


2 
—7 3 Esa bila T berada di antara P dan O 
2 2P0 


15 


3 5» 
— 0 bila PT kebalikan arahPO atau T disisi kiri P (5.109) 


(3,1) — Aa .In 3) (5.110) 
W4,1) —0 (5.111) 
(5,1) “ea (5.112) 
I(6,1) —0 (5.113) 
Untuk kasus dimana T berhimpit dengan O 
K1)--3.p04.r01mro- Iro| mro-13) (5.114) 
I(2,1)-0 (5.115) 
(3,1) — tak- terdefinisi (5.116) 
(4,1) — (5.117) 
(5,1) — tak- terdefinisi (5.118) 
I(6,1)- 0 (5.119) 
Untuk kasus dimana T berhimpit dengan P 
14,1) — 1 po 41 p0.mpo - 3ro(mro-1) (5.120) 
4 2 2 2 
12,1) -1.1.P0 (5.121) 
(3,1) — —1 (5.122) 
14,1)- 0 (5.123) 
5,1) — tak- terdefinisi (5.124) 
16,1) -— 0 (5.125) 
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Suatu kasus khusus lain yang penting untuk dipelajari adalah dimana titik T berada di 
tengah-tengah panel PO. Untuk kondisi tersebut rumus koefisien pengaruh dapat 
disederhanakan menjadi bentuk berikut 





I(41,0)--PO spon 2) -PO(inPO—1-In2) (5.126) 
(2,0) —7.7.PO (5.127) 
I(3,0)—0 (5.128) 
(4,0) —0 (5.129) 
4 
US, 0) 9 (5.130) 
I(6,0)—0 (5.131) 
1 1 PO 1 pa 
14,1) -— Ipa Hp. 2)- 1rofmro 1-In2) (5.132) 
(2,1) —7.x.PO (5.133) 
X3,1) ——1 (5.134) 
(4,1) -0 (5.135) 
KS) 5 (5.136) 
K6,1) H9 (5.137) 


Dari hasil-hasil yang diperoleh dalam pembahasan di atas kita dapat mengamati dan 
mengambil beberapa kesimpulan sebagai berikut : 

Untuk panel sumber kontinyu nilai koefisien pengaruh untuk fungsi potensial, fungsi arus 
dan komponen kecepatan v itu terdefinisi di semua titik di medan aliran tanpa 
perkecualian. Koefisien pengaruh untuk komponen kecepatan u juga mempunyai nilai di 
semua titik kecuali di titik-titik ujung panel dimana nilai koefisien pengaruh u menjadi tak 
terdefinisi. Koefisien pengaruh panel vorteks kontinyu pada umumnya juga terdefinisi 
kecuali untuk komponen kecepatan v di titik-titik ujung panel. Untuk panel doblet 
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kontinyu fungsi potensial dan komponen kecepatan u juga menjadi tak-terdefinisi di titik- 
titik ujung panel. 

Pengamatan tersebut di atas perlu diingat dengan sebaik-baiknya sebab dalam 
penyelesaian masalah dengan pendekatan metoda panel ada kemungkinan bahwa kita 
ingin menghitung nilai koefisien pengaruh di titik-titik ujung panel. Sekarang marilah kita 
periksa lebih teliti nilai koefisien pengaruh pada titik-titik di sepanjang garis yang tegak 
lurus (normal) pada panel PO dan melewati titik tengah panel. Untuk kasus ini rumus 
koefisien-koefisien pengaruh adalah 


K1,0) PO PO-In (7 PO TR? 4 TR.tan ag) 
X2,0) — ro| an) 4 aa 225) 
2 PO -PO 


X3,0)—0 
0..al. PO.TR 
(4,0) — tan en “2g a) 
4 
XS,0) 9 
N6,0)—0 


Variabel TR bernilai positif bila T berada di sebelah "atas" garis PO dan negatif bila 
berada dibawahnya. Bila titik T pada awalnya berada di atas PO dan kemudian 
diturunkan secara bertahap sampai berada di bawah PO maka nilai koefisien pengaruh 
pada umumnya juga berubah secara bertahap, Kecuali pada saat melintasi garis PO 


karena tanda TR berubah dari positif ke negatif. Suku-suku PO dan In (P0? #TR? 


tidak mengalami perubahan tanda bila TR berubah tanda. Tetapi perhatikan bahwa 





daa 2TR Etan 2TR 51 untuk TR positif 
P9 PO 





tan” 2TR #tan” Aa T untuk TR negatif 
PO -—PO 
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Ini berarti bahwa I(2,0) mengalami loncatan ("jump") sebesar t.PO bila titik T melintasi 
garis PO. Perhatikan pula bahwa sudut A, dimana 


A tan Psk ) 
TR? —PO?/4 


juga mengalami loncatan nilai sebesar 21 bila TR berubah dari nol positif (T sedikit di 
atas PO) menjadi nol negatif (titik T sedikit di bawah PO). Ini berarti bahwa nilai 1(4,0) 
mengalami loncatan sebesar 21 bila T melintasi garis PO. Nilai-nilai I(1,0), 13,0), 








(5,0), dan I(6,0) berubah secara mulus dengan perubahan nilai TR dan tidak mengalami 
loncatan pada saat T melintasi garis PO. Walaupun pembahasan di atas dibatasi pada 
kasus dimana T berada pada garis normal yang melewati titik tengah panel PO, namun 
kesimpulan yang bisa diambil juga berlaku untuk garis-garis normal lainnya. 

Sekarang marilah kita terapkan hasil pengamatan di atas untuk koefisien 
pengaruh fungsi potensial (lihat persamaan (5.14), (5.16) dan (5.18)). Fungsi potensial 
sumber, vorteks dan doblet masing-masing diwakili oleh koefisien I(1), X2) dan I(3). 
Koefisien I(1) dan I(3) bervariasi secara mulus sedangkan I(2) mengalami loncatan nilai 
bila melintasi garis PO. Ini berarti bahwa nilai fungsi potensial mengalami loncatan bila 
melintasi panel vorteks tetapi bervariasi secara mulus melintasi panel sumber atau panel 
doblet. 

Komponen kecepatan v diberikan oleh rumus-rumus (5.28), (5.30) dan (5.32). 
Nilai v untuk sumber, vorteks dan doblet masing-masing diwakili oleh koefisien I(4), 
X3), dan I(6). Pengamatan sebelumnya menunjukkan bahwa nilai I(4) mengalami 
loncatan bila melintasi panel PO sedangkan I(3) dan N(6) bervariasi secara mulus. Ini 
berarti bahwa kecepatan v untuk panel sumber mengalami loncatan bila melintasi panel 
PO. Secara fisis ini dapat dimengerti dengan mudah karena fluida memancar keluar ke 
arah atas (positif) di bagian atas PO dan memancar ke bawah (negatif) di bagian bawah 
PO. Nilai komponen kecepatan u untuk sumber, vorteks dan doblet diberikan oleh 
persamaan (5.27), (5.29) dan (5.31), yaitu masing-masing oleh W3), W4) dan XS). 
Pengamatan sebelumnya menunjukkan bahwa nilai u untuk panel vorteks mengalami 
loncatan bila melintasi panel PO sedangkan untuk sumber dan doblet nilai u bervariasi 
secara mulus. Variabel u dan v adalah komponen kecepatan tangensial dan normal pada 
panel PO dan loncatan nilai u melintasi panel PO menunjukkan bahwa sirkulasi pada 
sirkit tertutup yang mengurung elemen panjang panel vorteks PO tidak bernilai nol. Ini 
sesuai dengan pengertian kita bahwa y adalah kekuatan vorteks per satuan panjang dan 
sirkulasi itu sama dengan kekuatan vorteks. Oleh karena itu sirkulasi pada sirkit yang 
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mengurung elemen As dari panel vorteks PO adalah y dikalikan As. Kenyataan ini perlu 
diingat karena merupakan salah satu konsep dasar dalam penerapan metoda panel 
vorteks untuk menganalisa dan merancang bentuk aerofoil. Lembaran bidang datar 
singgularitas atau panel yang diwakili oleh garis lurus PO adalah balok-balok bangunan 
yang dapat dikombinasi secara linier untuk menyelesaikan masalah aliran potensial, tak 
mampat dengan syarat batas yang rumit. Sebelum kita membahas metoda panel dan 
penerapannya untuk analisis aerodinamika bentuk aerofoil, kita harus mempunyai 
pengertian yang mendalam mengenai sifat dan perilaku panel-panel sederhana yaitu panel 
tunggal baik untuk singgularitas sumber, vorteks ataupun doblet. Pola-pola aliran di 
sekeliling panel tunggal tersebut, baik di dalam aliran angin seragam ataupun di medan 
tanpa aliran akan dipelajari secara lebih terperinci dalam bab berikut. 
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LAMPIRAN 5.1 


Pengaruh panel singgularitas PO pada nilai fungsi potensial, fungsi arus dan aliran 
kecepatan di titik T tergantung pada koordinat titik-titik P, 9 dan T dan pada distribusi 
kekuatan singgularitas per satuan panjang. Nilai fungsi yang terimbas di T adalah hasil 
perkalian kekuatan dengan koefisien pengaruh panel PO. Nilai koefisien pengaruh 
ditentukan oleh geometri panel PO dan letak T. Bila koordinat P, O dan T masing- 
masing adalah (XP, YP), (XO, YO) dan (XT, YT) maka nilai koefisien pengaruh adalah 
fungsi dari beberapa parameter geometris yang dirumuskan di bawah ini. 

AXOP — XO- XP 
AYOP —- YO- YP 
AXTP —- XT-— XP 
AYTP — YT- YP 
AXTO —- XT-XO 
AYTO- YT-YO 


PO - (AXOP?  AYOP?)” 


PT — (AXTP? 4 AYTP)” 


OT- (AXTO? #AYTO?)” 
PR — (AXOP.AXTP 4 AYOP.AYTP)/PO 
-( 


TR — (AXOP.AYTP — AYOP.AXTP )/PO 
00- a28) 


TR 
00 — tan'| — 
O-tan (2) 





A-00 op tan PO-TE ) 


PT'—PO.PR 

Rumus-rumus koefisien pengaruh untuk kasus distribusi kekuatan linier diberikan 
di bawah ini, baik untuk fungsi potensial, fungsi arus maupun komponen kecepatan dan 
untuk singgularitas jenis sumber, vorteks serta doblet. 


Singgularitas Sumber/Serap 


Fungsi Potensial: — &(T,PO) -— Hreai.0) 4(70- "P).1(1,1) 
T 


Fungsiarus :  w(T,PO)- s-ler.X2,0) t(T0-—TP).1(2,1)| 
T 
Kecepatanu : u(T,PO) z 5-ler.X3,0) 4(70-— TP).1(3, 1)| 
Tt 
Kecepatan v : v(T,PO) - 5-ler.ita,0) 4(t0-1P).14, 
T 
Singgularitas Vorteks 


Fungsi Potensial: 6 T,PO)- — re, 0) #(y0 —yP).1(2, 1)| 


Fungsiarus :  w(T,PO)- — 1210) 410 -—P).KD)) 
T 
Kecepatanu :  w(T,PO)- — laa) #(y@—yP).1(4, 1)| 
T 
Kecepatanv :  v(T,PO)- sale. X(3,0) #(yO— YP). K3, 1) | 
T 
Singgularitas Doblet 


Fungsi Potensial: — &(T,PO)- sur. X(3,0) # (HO—uP).X(3,1)| 


Fungsi arus: w(T,PO) - —lnp.X(6,0) (HO - up). 14, 1)| 
Kecepatanu :  u(T,PO)- —lne.X5,0) #(nO-—pp).X(5,1)) 
Kecepatanv :  v(T,PO)- ——lnp.X(6,0) 4 (1O-uP). 16,1) 


Rumus umum integral I(n,k) 
T(4,1)-01(1,0) —-PO-OR.In OT-4-PR.In PT4TR.A 


X2,0) —-TR.In (OT/PT)- OR.00--PR.OP 
X3,0) —-In (OT/PT) 
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N4,0) — A 





























OR PR 
(5,0) — 2 
(5,0) OT PT 
TR TR 
X6,0) — . 
(6,0) OT? PT? 
Dn 2 2 7 Ag 2 
14,1) - PO 2PR | PO'-PR'4TR Ia OTP TR PT 4PR-TR 
4 2PO 2PO PO 
2. 2 2 Di 2 
K2,1) - IR PRIR (Oi ,PO-PR#TR 001PR TR' 99 
2 PO PT 2PO 2PO 
X3,1) IE Mn jeakLA 
PO LPT/ PO 
Kasi) Tam SJ iparA 
PO LPT/ Po 
Nag pa net, 
OT PO (PT 





TR A 
6,1) — — — 
Oa “0 
Integral I(n,k) untuk T disepanjang garis PO 
Y1,0) -—PO - OR.In|OR|- PR. In|PR| 


— 5» 
I(2,0) —0 bila OT searah dengan PO atau T disisi kanan O. 
T.OT bila T berada diantara P dan O 


T.PO bila PT kebalikan arah PO atau T disisi kiri P. 
X3,0) — —In((OR/PR)) 
I(4,0)—0 
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PO 
PR.OR 


I(6,0) —0 


(5,0) — 





PO -H-2PR 0 — PR? 
4 2PO 








K,1) — w|ORI:-. AP 


» D 
Y(2,1)-0 bila OT searah dengan P@ atau T disisi kanan O. 


K3,1) ——1-—In(|OR/PR)) 
I(4,1)—0 


1 1 

15,1) 5 —-—4-—I R/PR 

(5,1) -— Ga pom (OR/PR)) 
X6,1)-—0 

Integral I(n,k) untuk T berhimpit dengan O 


X1,0) — PO(In(PO) —1) 
— tak - terdefinisikan 


— tak - terdefinisikan 


) 

,1)— tak - terdefinisikan 
) 

,1)- tak - terdefinisikan 


(6,1) -—0 
Integral Y(n,k) untuk T berhimpit dengan P 
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Y(4,0) — PO(In(PO) —1) 
Y2,0) — xPO 
Y(3,0)— tak - terdefinisikan 
(4,0) —0 

(5,0) — tak - terdefinisikan 
(6,0) —0 


14,1) — PO(In(PO)-) 


I(6,1)-—0 
Integral I(n,k) untuk T di tengah-tengah panel PO 
I(1,0) - PO(In(PO) —1-In(2)) 


(2,0) — P0 
I(3,0) —0 
(4,0) —0 
(5,0) -—4/PO 
X6,0)—0 


14,1) -— 2PO(in(PO) -1-In(2)) 


(2,1) - srPO 
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BAB VI 
POLA ALIRAN PANEL SINGGULARITAS KONTINYU 


Dalam modul pertama bahan kursus ini kita telah membahas tentang pola aliran 
yang terimbas oleh singgularitas diskrit. Sifat linieritas persamaan Laplace dapat 
dimanfaatkan untuk memperoleh pola aliran rumit yang memenuhi kondisi batas 
tertentu, sebagai hasil kombinasi linier beberapa pola aliran elementer dari masing- 
masing singgularitas diskrit tersebut. Misalnya saja, pola aliran angin seragam yang 
menyapu bentuk lingkaran dapat diperoleh sebagai hasil kombinasi linier fungsi potensial 
kompleks untuk angin seragam dengan sebuah doblet dan vorteks diskrit yang berada 
atau terletak di titik pusat lingkaran. Pembahasan yang cukup rinci tentang pola aliran 
singgularitas diskrit telah diberikan dan sangat bermanfaat dalam penyelesaian masalah 
aliran fluida potensial di sekeliling bentuk lingkaran dan juga di sekeliling bentuk 
aerofoil. 

Panel singgularitas kontinyu dapat dianggap sebagai solusi elementer persamaan 
Laplace dan dapat digabungkan dalam suatu kombinasi linier untuk memperoleh solusi 
masalah aliran fluida potensial dengan kondisi batas yang rumit sesuai keinginan. Inilah 
sebabnya mengapa kita perlu membahas masalah pola aliran potensial yang diimbas oleh 
panel singgularitas kontinyu. 

Pola aliran yang diimbas oleh panel sumber, vorteks dan doblet akan dibahas 
secara berturut-turut, dimulai dengan kasus sederhana dimana kekuatan singgularitas per 
satuan panjang terdistribusi secara seragam dan dilanjutkan dengan kasus distribusi 
linier. 

Rumus-rumus untuk fungsi potensial, fungsi arus dan komponen kecepatan u 
dan v yang terimbas dititik T oleh panel sumber PO telah diturunkan dibab 5 dan ditulis 
ulang di bawah ini. 


dx) .Flx,y) (6.1) 

T 

vtxsy) ——— .F2lx,y) (6.2) 
IT 

u(x,y) .F3(x,y) (6.3) 
T 


v(x,y) - Fata) (6.4) 


F1(x,y) —-PO- OR (x,y).In OT(x, y) “PR (x,y)-In PT(x,y) #TR(x, y). A(x,y) (6.5) 








F2(x,y) — “ebay M3) — OR(x,y).@0(x,y) #PR (x,y).OP(x,y)(6.6) 
PT(x,y) 

F3(x,y) - al Ona) ) (67) 

F4(x,y) ad A(x,y) (6.8) 


Kuantitas & adalah konstanta yang mewakili kekuatan sumber per satuan panjang, 
sedangkan PO, PR, TR dan seterusnya dirumuskan sebagai berikut 

















po- ((xo-xp) 4(y0-—ve))” (panjang panel) (69) 
PR(x,y) - ((XO-XP).(x— XP) 4 (YO- YP).(y— YP))/PO (6.10) 
OR(x,y) -PR(x,y)-PO - ((XO- XP).(x— XP) (YO-— YP).(y— YP))/PO 

(6.11) 
TR(x,y)- ((XO-XP).(y-YP)-(YO- YP).(y- YP))/P9 (6.12) 
PT(x,y) - ((x- xp) 4 (y— xe))? (6.13) 
@rts,y) - (x-xo) 4(y- vo))” (6.14) 
00(x,y) - an sat (6.15) 
OP(x,y) — aa) (6.16) 

ba PO.TR (x,y) 

Alx,y) -00(x,y) -6P(x,y) — tan bara 4 (6.17) 


dimana (XP,YP), (XO,YO) dan (x,y) masing-masing adalah koordinat dari titik P, O, 
dan T. 
Perhatikan bahwa 
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FLGy) - “tny) (618) 


jadi F1(x,y) sesungguhnya adalah fungsi potensial non-dimensional. Ini berarti bahwa 
kurva isopotensial adalah letak atau kedudukan dari titik-titik dimedan aliran 
disepanjang mana fungsi F1(x,y) bernilai tetap. Koordinat (XK, YK) dari titik-titik kurva 
isopotensial dengan nilai fungsi potensial non-dimensional K adalah solusi dari 
persamaan f(x,y)-0 dimana 

f(x,y)- Filx,y) —K (6-19) 


Untuk suatu nilai x tertentu, misalnya x-XK, persamaan (6.19) dapat disederhanakan 


menjadi 

G(y) -f$ (XK,y) - FIXK,y) -K (6.20) 
dan YK adalah solusi dari persamaan 

G(y)-0 (6.21) 
Sebaliknya apabila nilai y-YK yang ditentukan maka persamaan (6.19) dapat ditulis 
menjadi 

H(y) f(x, YK) —- Fix, YK)-K (6.22) 
dan XK adalah solusi dari persamaan 

H(y)-0 (6.23) 


Persamaan (6.21) dan (6.23) adalah persamaan aljabar dalam satu variabel bebas, yaitu 
variabel real x atau y, yang solusinya dapat dicari dengan menggunakan salah satu dari 
metoda-metoda numerik seperti metoda pembelah selang ("bisection"), metoda posisi 
palsu ("regula falsi"), metoda sekan, metoda Newton-Raphson dan lain sebagainya. 

Prosedur komputasi koordinat titik-titik kurva isopotensial dapat dijelaskan 
sebagai berikut. 

Pertama-tama batasan medan solusi harus ditentukan, misalnya berbentuk 
segipanjang ABCD dimana koordinat titik A, B, C dan D adalah (XMIN,YMIN), 
(XMAK,YMIN), (XMAK,YMAK) dan (XMIN,YMAK). Suatu jaringan kisi 
komputasional yang bersifat ortogonal (saling tegak lurus) dapat dibentuk dengan 
"menggambar" famili garis vertikal berindeks i dan famili garis horisontal berindeks j 
(lihat gambar (6.1)). Nilai indeks i ditentukan bervariasi dari 1 sampai dengan JMAX. 
Rumus untuk garis kisi vertikal berindeks i adalah 


x(i) —- XMIN --(i—1).Ax untuk 1 Si K IMAX (6.24) 
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sedangkan rumus untuk garis kisi horisontal adalah 

y(i) — XMIN #(j —1).Ay untuk 1 « j « JMAX (6.25) 
Selang Ax dan selang Ay biasanya ditentukan bernilai sama, misalnya h. Nilai XMIN, 
YMIN, h, IMAX dan JMAX diberikan sebagai data masukan ("input"), sedangkan nilai 
XMAX dan YMAX dihitung sebagai berikut 

XMAX — x(IMAX) — XMIN 4-(IMAX -1).h (6.26) 

YMAX — y(JMAX) — YMIN -(JMAX-1).h (6.27) 
Perpotongan antara garis x(i) dengan garis y(j) disebut titik kisi T(i,j) dengan koordinat 
(x@),yG)). Karena koordinat T(i,j) diketahui maka nilai fungsi F1(x,y) di titik T(i,j) dapat 
dihitung dengan rumus (6.5) dengan hasil misalnya F1(,j). Untuk suatu nilai fungsi 
isopotensial tertentu, misalkan K, kita dapat mencari solusi dari persamaan (6.19), yaitu 
fx.y) — 0, berupa koordinat titik yang berda pada garis kisi vertikal x(i) dengan 


menyelesaikan persamaan 

Gli,y) -#(xi),y) —0 (6.28) 
Solusi berupa titik pada garis kisi horisontal y(j) dapat dihitung dengan menyelesaikan 
persamaan 


H(jx) f(x y(i)) 0 (6.29) 
Nilai K tidak boleh ditentukan sesuka hati dan harus memenuhi persyaratan tertentu. 
Untuk masing-masing garis kisi vertikal x(i) nilai-nilai fungsi F1(i,j) harus disortir untuk 
menentukan nilai minimal dan maksimalnya, misalnya FYMIN(i) dan FYMAXII), 
sedemikian rupa sehingga Untuk setiap nilai i dimana 1 Xi « IMAX 

FYMIN(i) SFXi, j) SFXMAX(i) untuk 1 « j « JMAX (6.30) 
Selanjutnya nilai FYMIN(i) harus disortir untuk menentukan nilainya yang terbesar 
misalnya FYMINMAX, dimana 

FYMIN(i) KS FYMINMAX untuk 1Sis IMAX (6.31) 
Nilai-nilai FYMAXG) juga harus disortir untuk mencari nilainya yang terkecil, misalnya 
FYMAXMIN, dimana 

FYMAX(i) 2 FYMAXMIN untuk 1 Si S IMAX (6.32) 
Persamaan (6.28) hanya mungkin mempunyai solusi pada garis kisi x(i) apabila K 
memenuhi persyaratan berikut 

FYMIN(i) KK SFYMAXII) 

(6.33) 
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Untuk masing-masing garis kisi horisontal y(j) nilai-nilai fungsi F1(i,j)) harus disortir 
untuk menentukan nilai minimal dan nilai maksimalnya, misalnya FXMING) dan 
FXMAX(), sedemikian rupa sehingga untuk setiap nilai j dimana 1 £ j X JMAX 
FYMINMAX £ K « FYMAXMIN (6.34) 

Selanjutnya nilai-nilai FKMING) dan FXMAX()), harus disortir untuk menentukan nilai- 
nilai FKMINMAX dan FMAXMIN yang memenuhi persyaratan berikut 

FXMIN(i) S Fl(i,j) SEKMAX untuk 1 Si £ IMAX (6.35) 
Solusi dari persamaan (6.29) hanya mungkin ada pada garis horisontal y(j) apabila K 
memenuhi persyaratan 

FXMIN(j) SEXMINMAX untuk 1 j  JMAX (6.36) 

FXMAX(j) 2 FKMAXMIN untuk 1£ js JMAX (6.37) 
Solusi dari persamaan (6.29) hanya mungkin ada pada garis horisontal y(j) apabila K 
memenuhi persyaratan 

FXMIN(j) SK S FXMAX/j) (6.38) 
Pertimbangan untuk seluruh medan solusi menunjukkan bahwa sebaiknya nilai K 
memenuhi persyaratan berikut. 

FXMINMAX £ K « FXMAXMIN (6.39) 
Perlu diperhatikan bahwa persyaratan (6.34) dan (6.39) telah diperoleh berdasarkan 
suatu pertimbangan tertentu yang mungkin tidak sesuai distribusi nilai-nilai FIi,j) 
dimedan solusi. Ada kemungkinan bahwa distribusi tersebut bersifat sedemikian rupa 
sehingga 

FXMINMAX » FXMAXMIN 
Sehingga persyaratan (6.34) menjadi mustahil untuk dipenuhi. Suatu persyaratan untuk 
nilai K yang lebih longgar adalah dengan mencari nilai terkecil, FIMIN, dan nilai 
terbesar, FIMAX, dari fungsi F1(i,j)) untuk semua nilai i dan j diseluruh medan solusi 
untuk menuntut agar 

FXMINMAX » FXMAXMIN 
Walaupun persyaratan (6.40) dipenuhi, ada kemungkinan bahwa hanya sebagian dari 
kurva isopotensial yang berada dimedan solusi sedangkan bagian lainnya berada diluar 
medan solusi sehingga kita tidak dapat "menggambar" kurva isopotensial seutuhnya. 

Perlu diingat bahwa bentuk kurva yang ingin "digambar" itu tidak diketahui 
sebelumnya dan dapat berwujud garis yang pada dasarnya horisontal dari kiri ke kanan, 
atau garis yang pada dasarnya vertikal dari atas ke bawah., dan mungkin juga berbentuk 
suatu kurva tertutup yang rumit (lihat gambar (6.2), (6.3) dan (6.4)) 
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Setelah batasan medan solusi dan nilai K ditentukan, kemudian kita dapat 
menghitung koordinat titik-titik pada kurva isopotensial sebagai berikut. Untuk setiap 
nilai i, dimana 1XiX IMAX, kita cari solusi dari persamaan (6.28) sebagai berikut. 
Nilai y dimana G(i,y)-0 pasti berada diselang y(j) Sy S y(j #1), apabila persyaratan 
berikut dipenuhi 

GM, )-.GGj--1) co (6-41) 
Oleh karena itu untuk setiap nilai j dimulai dari j-1 sampai j-JMAX-1 nilai Gi,j). 
Gfi,j#l) dihitung dan apabila persyaratan (6.41) dipenuhi misalnya untuk j-JS maka 
diketahui bahwa ada solusi digaris x(i) dan solusi tersebut berada diselang 
Y(IS)YSY(JS-H1) dimana YS adalah solusi yang dicari. 

Setelah nilai i ditambah 1 dan prosedur diulang lagi untuk mencari solusi garis 
x(i) berikutya. Apabila kurva yang ingin digambar pada dasarnya berbentuk kurva 
horisontal maka disetiap garis x(i) pasti ada solusi yaitu koordinat (XS(i),YS(i)) dimana 


XS(i) — xi) (6.42) 
dan YS(i) adalah solusi dari persamaan 
G(i,y)-0 (6.43) 


Apabila kurva yang ingin digambar pada dasarnya berbentuk kurva vertikal maka 
prosedur penyelesaian masalah harus disesuaikan yaitu diubah sebagai berikut. Untuk 
setiap garis kisi horisontal atau untuk setiap nilai j, dimana 1X jX JMAX, pasti ada 1 
solusi yaitu titik dengan koordinat (XS), YSG)) dimana 


YS(H) -yY(H) (6.44) 
dan XS(j) adalah solusi dari persamaan 
H(j,x) —0 (6.45) 


Letak selang i dimana solusi XS berada dapat diperoleh dengan menerapkan persyaratan 
bahwa diselang tersebut 

H(i))-H@ 4, j) «0 (6.46) 
Sekarang perhatikan bahwa prosedur diatas tidak dapat diterapkan untuk kasus yang 
bersifat lebih umum dimana pada setiap garis vertikal x(i) atau garis horisontal yg) 
belum tentu selalu ada 1 solusi, tetapi mungkin tidak ada solusi atau ternyata ada banyak 
solusi. 

Untuk kasus umum tersebut kita harus menggambar kurva yang diinginkan 
sepotong-sepotong, di setiap elemen segi-panjang kisi yang dibatasi oleh titik-titik T(i,j), 
Tj), TGtLj 1) dan TM,j#1) untuk 1KIKIMAX-—1 dan IX jX JMAX-1. 
Untuk setiap nilai i dan j, segi-panjang Tlij), Wi#1,)), TG#1,j H1) dan T(i,j-1) disebt 
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sebagai kotak (i,j) dengan empat sisi yaitu TG )TGj#1), T1)T(L,j-1), 
TG-#Lj#-)TG-#L,j) dan T(itL,j)Tfi,j) sebagai sisi pertama, kedua, ketiga dan keempat 
kotak tersebut. Selanjutnya kita membuat anggapan bahwa pada kotak (i,j) ada 
kemungkinan bahwa 2 dari keempat sisi mempunyai titik solusi. Tetapi ada juga 
kemungkinan tidak ada solusi pada ke 4 sisi kotak tersebut. 
Logika prosedur hitungan demikian dapat dijelaskan sebagai berikut. 
1. Rumuskan nilai awal JS-0. 
2. Uji apakah ada solusi di sisi pertama, yaitu apabila 
KG)aGj td) co 
maka ada solusi di sisi pertama. Kemudian nilai JS ditambah 1, yaitu 
JIS—JS-1-1 
dan koordinat solusi dapat dihitung dan dicatat sebagai (XS(JS), YS(JS)) dimana 
XS(JS) — xli) 
sedangkan YS(JS) adalah solusi dari 
G(i,y)-0 
apabila ffi,j).fi5j#1)-0 maka f(i,j)-0 atau f(i,jt1)—0. Untuk kasus f(i,j))-0 dan 
f(i,j #1) & 0 maka 


XS(JS) — xi) 
YS(JS) — y(i) 
sedangkan bila f (i,j) #0 dan f(i,j-1)-0 maka 


XS(JS) — x(i) 

YS(JS) — y(j #1) 
Sebaliknya apabila ffi,j).i,j-1)20 maka tidak ada solusi pada sisi pertama kotak 
(i, j) dan JS tetap bernilai 0. 

3. Uji apakah ada solusi disisi kedua. Apabila fti,j#1)-0 maka titik T(i,j#1) adalah 
solusi, tetapi sudah tercakup dalam butir 2 sebagai solusi disisi pertama dan tidak 
boleh diangap solusi lagi untuk kedua kalinya. Ini berarti bahwa jika (i,j-1)-0 
maka tidak ada solusi disisi kedua. Apabila f(i#1,j#1)-0 maka titik TG#1,j#1) 
adalah titik solusi disisi kedua. Demikian juga apabila fi,j).fi,j#1)XO maka ada 
titik solusi pada selang horisontal antara T(i,j#1) dan T(iH1,jt1) yang nilainya 
dapat dicari dengan menyelesaikan persamaan H(j-1,x)-—0. 

Untuk kasus dimana ada solusi disisi kedua, nilai JS harus ditambah 1 sehingga 


menjadi 
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JS—1 bila dibutir 2 JS bernilai 0 
atau JS-2 bila dibutir 2 JS bernilai 1 

Apabila nilai baru JS adalah 2 maka garis solusi yang menghubungkan 

(XS(1),YS(1)) dan (XS(1),YS(2)) harus digambar dan kita kembali ke butir 1 

untuk menghitung solusi dikotak (i,j) berikutnya. Bila JSX2 kita teruskan 

prosedur ke butir 4. 

4. Uji apakah ada solusi di sisi ketiga. 

Apabila ffit1,j#1)-0 maka di sisi ketiga tidak ada solusi karena T(i#1,jH1) 

dianggap sebagai solusi di sisi kedua. Apabila fti#1,j)/-0 maka titik T(i#1,j) 

adalah titik solusi di sisi ketiga. Demikian juga bila fi--1,j#1).fi#1,j)X0O maka 

ada solusi pada selang vertikal antara T(i-1,j-1) dan T(i#1,j) yang nilainya dapat 

dicari dengan menyelesaikan persamaan G(i-1,y)-fi#1,j-#-1).f--1,j)-0. Untuk 

kasus dimana ada solusi di sisi ketiga, nilai JS harus ditambah 1 sehingga menjadi 
JS — 1 bila dibutir 3 JS bernilai 0 

atau JS —2 bila dibutir 3 JS bernilai 1 

Apabila nilai baru JS adalah 2 maka garis solusi yang menghubungkan 

(XS(1),YS()) dan (XS(2),YS(2)) kemudian digambar dan kita kembali ke butir 

1 untuk menghitug solusi di kotak (i,j) berikutnya. Jika JS#2 maka kita 

langsung ke butir 1. 

Dalam prosedur diatas telah dijelaskan bahwa YS(JS) adalah solusi dari 
persamaan G(I,y)-0 atau XS(JS) adalah solusi dari persamaan H(J,x)-0 dan lain 
sebagainya. Solusi dari persamaan tersebut dapat diperoleh dengan menggunakan 
metoda numerik seperti metoda Newton. walaupun metoda Newton itu cukup efisien 
tetapi upaya komputasi yang dibutuhkan masih tetap besar, karena fungsi G(y) atau H(x) 
dan turunannya, yaitu G'(y) atau H'(x), berbentuk cukup rumit. 

Suatu pendekatan lain yang dapat diterapkan adalah dengan menaksir nilai solusi 
dengan metoda interpolasi. Metoda ini sangat efisien walaupun keakuratannya terbatas. 
walaupun demikian keakuratan hasil taksiran dapat dianggap cukup baik, khususnya bila 
rumus interpolasi yang digunakan adalah polinom derajat 2 atau derajat lebih tinggi. 

Prosedur yang diberikan di atas memang telah diberikan untuk penggambaran 
kurva isopotensial pada khususnya, tetapi sesungguhnya berlaku umum untuk 
menggambar kontur kurva iso-nilai fungsi apapun juga. 

Selanjutnya perlu diperhatikan bahwa fungsi F1(x,y) yang diberikan oleh rumus 
(6.5) sesungguhnya tidak hanya mewakili fungsi arus singgularitas vorteks. 
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Fungsi F2(x,y) yang diberikan oleh persamaan (6.6) disamping mewakili fungsi 
arus untuk sumber juga mewakili fungsi potensial untuk vorteks. F3(x,y) yang diberikan 
oleh persamaan (6.7) mewakili fungsi potensial untuk doblet disamping. mewakili 
komponen kecepatan u untuk sumber dan komponen kecepatan v untuk vorteks. 
Selanjutnya fungsi F4(x,y), persamaan (6.8), mewakili fungsi arus untuk doblet, 
komponen kecepatan v untuk sumber dan komponen kecepatan u untuk vorteks. 

Komponen kecepatan u untuk doblet diberikan oleh rumus berikut 


OR(x,y) PR(x,y) 





F5(x,y) — 6.47 

ON Lanny PTN in 
sedangkan kecepatan v untuk doblet diberikan oleh rumus 

F6(x,y) — TR(x,y) 3 TR(x,y) (6.48) 


OT(x,y) PT(x,y) 
Dengan pengertian di atas kita dapat mempelajari kontur-kontur isopotensial, garis arus, 
kecepatan u dan kecepatan v untuk semua singgularutas, yaitu sumber, vorteks dan 
doblet, hanya dengan mempelajari kurva-kurya iso-nilai fungsi F1, F2, F3, F4, F5 dan F6 
saja. 

Untuk kasus dimana kekuatan singgularitas diberikan oleh suatu distribusi linier, 
kita perlu mempelajari kurva-kurva iso-nilai fungsi-fungsi berikut. 


PO-4-2PR(x,y) PO'-—PR(x,y) #TR(x,y) 
3 y 























FIXx,y) — : AT In OT(x,y) 
, 2 (6-49) 

Pa 2 R3) MPT,y)4 PR IR) Sales 

Kota Tan) PRE ARA wa ten) 
(x.y) 
HER #TRY O0(x,y) ai ka Intan ny 

(6.50) 
F13(x,y) 5-1 2 y) Ja Iran). na .Nx,y) (6.51) 
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F14(x,y) Ba .i Intan ur 2 Nx,y) (6.52) 

F15(x,y) — @R(x,y) , . ES (6.53) 
OTx,y) PO PT(x,y) 

Fletx,y) 10), IR(R,y) (6.54) 


Po OTtx,y) 

Bentuk kurva iso-nilai fungsi F1(x,y), F2(x,y) dan seterusnya dapat dilihat pada gambar 
(6.5) sampai dengan (6.16). Karena persamaan Laplace bersifat linier maka kita dapat 
memperoleh solusi yang lebih rumit sebagai hasil gabungan linier 2 atau lebih fungsi- 
fungsi elementer. Dalam masalah aerodinamika pada dasarnya yang ingin dipelajari 
adalah aliran angin seragam yang diganggu oleh keberadaan suatu bentuk seperti 
lingkaran, aerofoil dan lain sebagainya. 

Masalah tersebut dapat dimodelkan sebagai aliran angin seragam yang diganggu 
oleh lembaran singgularitas. Untuk memperoleh gambaran tentang pengaruh keberadaan 
lembaran singgularitas pada pola aliran angin seragam, kita perlu mempelajari pola-pola 
aliran yang diperoleh sebagai hasil kombinasi linier antara angin seragam dan satu panel 
singgularitas. 

Fungsi potensial kompleks untuk angin seragam adalah 


O(x,y) - (x-tiy).(cosa—isina) (6.55) 
sedangkan turunannya adalah 

D'(x,y) —u-—iv — cosa—isina (6.56) 
Ini berarti bahwa fungsi isopotensial dan fungsi arus untuk angin seragam adalah 

6(x,y) —xcosa #ysina (6.57) 

w(x,y) — ycosa —xsina (6.58) 
sedangkan rumus untuk komponen kecepatan adalah 

u(x,y) 5 c9osa (6.59) 

v(x,y) —sina (6.70) 


Fungsi potensial angin seragam yang diganggu oleh sebuah panel sumber PO, diberikan 


oleh rumus berikut 
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0xy) -xeosat ysina 4-—-FXx,y) (6.71) 
T 


Fungsi arus untuk aliran yang sama adalah 








wl(x,y) — ycosa xsina t-—-FZx,y) (6.72) 
T 


sedangkan rumus untuk komponen kecepatan adalah 


u(x,y) — coso.t-—K3(x,y) (6.73) 
T 


v(x,y) -sina 4” Fdlx,y) (6.74) 
T 


Rumus-rumus diatas berlaku untuk panel sumber dimana kekuatan sumber terdistribusi 


secara seragam. 


Pola aliran yang diperoleh jelas tergantung pada nilai 22 Untuk kasus dimana aa «cl 
T T 


atau O x« 21 jelas bahwa aliran angin seragam hampir tidak berubah sama sekali 
sedangkan untuk & £ 21 pengaruh keberadaan panel sumber tentu sangat mengganggu 
pola aliran angin seragam. 

Bentuk pola aliran atau garis-garis arus untuk kasus dimana panel PO berada 
sejajar angin arah angin, dengan koordinat titik P dan 9 bernilai (-1/2,0) dan (1/2,0) 
dan &—0.21x ditunjukkan pada gambar (6.17), sedangkan untuk & — 2x ditunjukkan 
pada gambar (6.18). 

Pola aliran untuk kasus dimana panel sumber PO berdiri tegak lurus arah angin 
seragam ditunjukkan pada gambar (6.19) dan (6.20) dimana koordinat P dan 9 adalah 
(0,—1/2) dan (0,1/2). Seperti dapat dilihat, berbagai berbagai macam pola aliran dapat 
diperoleh sebagai hasil kombinasi linier angin seragam dengan panel sumber distribusi 
kekuatan seragam, hanya dengan mengubah kekuatan sumber dan orientasi atau sudut 
yang dibentuk oleh arah angin dan panel PO. Pola aliran yang lebih bervariasi dapat 
diperoleh dengan menggunakan distribusi kekuatan linier. Kekuatan per satuan panjang 
di titik P dan O ditentukan bernilaa oP dan 0G, sehingga distribusi kekuatan 


disepanjang panel PO diberikan oleh rumus 
als) -GP “(00 —0P).(s/PO) 
(6.75) 
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atau 
als) -GP(1—(s/PO)) -60.(s/PO) (6.76) 


Dari hasil pembahasan di bab 5 kita dapat menurunkan rumus-rumus untuk fungsi 
isopotensial, fungsi arus dan komponen kecepatan untuk kasus distribusi linier (6.76) 


sebagai berikut. 

@xy) - SP (ri(x,y) -F11x,y)) 4 Orii(x,y) (6.77) 

2x 21x 
wx,y) - (r2x,y) -FI2x,y)) #2 Or12(x,y) (6.78) 

2x 2x 

utny) SP (F3(xy) -F13(x,y)) SO r13(x,y) (6:79) 

21x 21x 
vaoy) «SP (rato,y) -F1A(x,y)) t Oria(x,y) (6.80) 

21x 2x 


Untuk kasus aliran angin seragam yang diganggu oleh panel sumber PO dengan 
distribusi kekuatan linier, rumus-rumus untuk pola alirannya adalah 


dlx,y) —x.c0SO #y.sino 1 rita) — Fii(x, y)) 1 ri Ux,y) (6.81) 
T Tt 








w(x,y) —y.cosa —x.sina 1 (ray) -Fiz(x,y)) 120 126,682) 

u(x, y) 5 cosa- og (F3(x, y) F13(x,y)) 100 p13(x,y) (6.83) 
2x 2m 

v(x,y) Sena (Fa(x,y) F1dtx,y)) - Oriax,y) (6.84) 
21 2x 


Seperti dapat dilihat persamaan-persamaan (6.81) sampai dengan (6.84) itu lebih rumit 
dibandingkan dengan (6.71) sampai dengan (6.74). Rumus untuk distribusi linier 
melibatkan satu konstanta tambahan dan oleh karena itu variasi pola aliran yang dapat 
diperoleh menjadi lebih banyak. 

Pola garis arus untuk kasus dimana oP-x1 dan 09 -——1, oP-—0.1x1 dan 
600——21x dan cP——2x dan 60 —0.11 ditunjukkan pada gambar (6.21), (6.22) dan 
(6.23) untuk kasus dimana koordinat tiik P dan 9 adalah (—1/2,0) dan (1/2,0). Untuk 


nilai-nilai cP dan 60@ seperti diatas tetapi dengan koordinat titik P dan O diberikan 
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bernilai (0,—1/2) dan (0,1/2), yaitu PO tegak lurus arah angin, bentuk pola aliran 
diberikan pada gambar (6.24), (6.25) dan (6.26). Rumus-rumus untuk aliran angin 
seragam yang diganggu oleh vorteks PO adalah sebagai berikut. 


blx, y) —x.c0SA #y.sin ad — 1 (rata, y) -F12(x,y)) 1 125,563) 





v(xy) —y.cosa—x.sina (ri6.y) F11(x,y))-— pit) 86) 
TI 





ulx,y) -cosa— YP. (Falx,y) —Fta(x,y)) —"Oriatx,y) (6.87) 
2x 2x 


v(x,y) — sin ad " (r3(x,y) F13(x,y)) 1S a33) (6.88) 





Pola garis arus untuk aliran ini ditunjukkan pada gambar (6.27), (6.28), (6.29), (6.30), 
(6.31) dan (6.32). 

Rumus-rumus untuk aliran angin seragam yang diganggu oleh panel doblet PO 
adalah sebagai berikut 


6(x,y) —x.coso #y.sina #E (ratan) -F13(x,y)) ita, y)(6.89) 





w(x,y) —y.cosa xsino— BE (rata) F1(x,y)) BE L4.3)630) 

akan) masa ena ah F15(x,y))- EO 15(x,y) (6.91) 
2x 2x 

v(x,y) — sin —"F (Felx,y) -F16(x,y)) 0 r1elx,y) (6.92) 
2x 2x 


Pola garis arus yang diperoleh dari rumus (6.90) untuk beberapa kombinasi linier 
UP, yO dan koordinat P dan 9 ditunjukkan pada gambar (6.33), (6.34), (6.35), (6.36), 
(6.37) dan (6.38). 

Seperti dapat dilihat berbagai macam pola aliran dapat diperoleh sebagai hasil 
kombinasi linier angin seragam dan satu panel singgularitas. 

Pola-pola aliran yang jauh lebih rumit tentu saja dapat diperoleh dengan 
memperbanyak jumlah panel singgularitas. Sebagai contoh marilah kita pelajari pola 
aliran angin seragam yang diganggu oleh empat panel vorteks, yaitu panel POI, PO2, 
PO3 dan PO4. Untuk memudahkan pembahasan ditentukan bahwa panel PO1 adalah 
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panel dengan titik-titik ujung P(1) dan O(1), yang berkoordinat (XP(1), YP(1)) dan 
(XO(1), YO(1)). Secara umum ditentukan bahwa panel ke J adalah panel PO(J) dengan 
titik-titik ujung P(J) dan O(J). Koordinat kedua titik tersebut adalah (XP(J), YP(J)) dan 
(XO(I), YOM)). Selanjutnya koefisien pengaruh fungsi arus untuk panel ke J pada titik 
T diberi simbol FI(T,J) dan F11(T,J), dimana 


FLT, J) — -PO(I)- OR(T, J).In OT(T, J) -PR(T,J)-In PTCT,J) 








(6.93) 
#TR(T,J).ACT,J) 
FIT, J) — H0) #2PR(T, J)) “ram BOHAI Ja OT(T, J) 
h HP RI Jpn, pa PRT)-TR(T,I) MTI) 
2 PO(J) POJ) 
(6.94) 
p0(3) - ((xo(3)- xe(3))” #vota)- ve) (6.95) 
Patra) KO) xp). (xT— xp) 1 (xo(a) - 3 MU AR 
(6.96) 
OR(T,J)-PO(T,J)-PO(J) (6.97) 
PT(T,3) Ma aa UR 
(6.98) 
pT(3) - ((xr—xp/(3)) : vr ve)? (6.99) 
013) - ((x1-x0(3)) Hr yo) (6.100) 
00(T, J) — Pat) (6.101) 
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OP(T, J) -— wa at) (6.102) 


A(T,J) -60(T,J) -6P(T, J) (6.103) 





(6.104) 





2 POJ).TR(T,J) 
A(T,J) — tan Ec Na er 


Rumus untuk fungsi arus dititik T adalah 


xp(4)—1/2/2 — YP(4)-1/2/2 — yP(4)—-1/2 
x0(4)-1/2/2  YO(4)--1/2V2 yO(4)-1/2 





W(T) - YT.cosa — XT.sina 2 2 eta) rta) yA) 


Zian pa AW LixT, n| 


Jumlah panel tidak perlu dibatasi bernilai 4 tetapi bernilai berapapun juga. Bila jumlah 
panel adalah JP maka rumus di atas dapat ditulis menjadi 


YT)- WT,as) HI Sora) BJ cra, J. rota)| (6.105) 


Te 
dimana 
W/(T) adalah nilai fungsi arus dititik T 
W(T,as) — YT.cosa — XT.sina adalah nilai fungsi arus dititik T karena berada dalam 
medan aliran angin seragam yang membentuk sudut « dengan sumbu horisontal 
CFP(T,J) - FXT, J)-F1LT, J) 
CFO(T,J) - FIXT,J) 
Gambar (6.39) menunjukkan pola aliran yang diberikan oleh persamaan (6.105) untuk 
data-data berikut 


xP(1)-1/2/2 — YP(1)--—1/2/2 yP(1)-1/2 
x0(1)--1/2/2 YO(1)--1/2/2 yO(1)- 
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xP(2)——1/2/2  Yp(2)--—1/2/2 yP(2)-1 
xX0(2)--1/2/2 “vO(2)-1/2/2  y0(2)--1 


xXP(3)——1/2/2  YP(3)-1/2/2  yP(3)—-1 
x0(3)-1/2/2 — v@O(3)-1/2/2 y0(3)--1/2 


xP/(4)-—1/2/2 — yP(4)-1/2/2 — yP(4)——1/2 
x0(4)-1/2/2  vOl4)--—1/2/2 y0(4)-1/2 


Gambar (6.40) menunjukkan pola aliran yang diberikan oleh persamaan (6.105) untuk 
data-data berikut 


xXP(1)-1  yP(1)-0  yP(1)-0 
x0(1)-0 yY@1)--4 y@(1)-2 
XP(2)-0 YP(2)-—1 ”yP(2)-2 
X0(2)-—1 YO(2)-0 y0(2)-0 


XP(3)-—1 YP(3)-0 yP(3)-0 
x0(3)-0  YO(3)-1 yO(3)-—2 


XP(4)-0  YP(4)-1 yP(4)-—2 

x0(4)-1 YOl4)-0 yOla)-0 
Dalam contoh diatas dapat diperhatikan bahwa titik O(1) berhimpit dengan P(2), 0(2) 
berhimpit dengan P(3), O(3) berhimpit dengan P(4), O(4) berhimpit dengan P(1). Secara 
umum dapat dikatakan bahwa titik O(J) berhimpit dengan P(JH1) untuk 1S JS JP-—1 
dan O(JP) berhimpit dengan P(1) sedemikian rupa sehingga kurva yang dibentuk oleh 
panel-panel tersebut merupakan suatu sirkit atau suatu kurva tertutup. 


Untuk kasus seperti ini persamaan (6.105) dapat ditulis ulang menjadi bentuk 
berikut 


YT)- WT,as)-— Hr CFP(T,J).yP(J) P3 crotr,3)ata)| 


J-1 J-1 


atau 
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T| Ja J-2 


YT)- WT,as)-— HS CFP(T,J).yP/(J) 1 crofr, J5 Data) 


WT)- WT,as) - jar, 1).yp(1) Pa (CFP(T, J).CFO(T,J —1)).y2/(J) 


J-2 
#CFO(T, JP). yB(JP 41)| 
Mengingat bahwa titik J-1 dan J-JP-#1 itu berhimpit maka 
YP(JP #1) —7P(1) 
dan persamaan diatas dapat disederhanakan menjadi 


VP) Was) Ca). rRt) (6.106) 
dimana 


C(T, 1) — CFP(T,1)-CFO(T, JP) 
5 FUT,1)-F1NT,1) #F1UT, JP) 
C(T,J) —- CFP(T,J)-CFO(T,J —1) untuk 2S J X JP 
5FKT,J)-F1UT,J) #F1XT, J -1) 
Untuk kasus dimana titik J-1 dan J-JP4#1 tidak berhimpit, rumus (6.106) diubah sedikit 
menjadi 
1 JP 
WT) - w(T,as) 2g 2 CT,3)1PG) (6.107) 
dimana 
C(T, 1) — CFP(T,1) - FXT,1)-F1NT,1) 
C(T,J) - FUT,J)-F1KT,J)#FIKT,J—1) untuk 2s Jx JP 
C(T, JP #1) — CFO(T, JP) — FIT, JP) 
Perhatikan bahwa dalam kita boleh menentukan nilai YP(J) dengan sesuka hati, tetapi 
kita tidak tahu bentuk pola aliran yang akan diperoleh dengan pilihan nilai-nilai YP(J) 
tersebut. 
Seandainya kita menganggap bahwa yP(J) adalah nilai-nilai yang tidak diketahui 
dan sebaliknya kita menganggap bahwa nilai w(T) disejumlah titik T(I), misalnya untuk 
1X IX JP, itu diketahui maka persamaan (6.106) dapat ditulis sebagai berikut. 
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WI) WTGjas) -——Y CM). rp) 


yang dapat ditulis ulang menjadi 


SH 04,3).yeta) - 21. W(T(I),as) —21. y(I) (6.108) 


J-1 
untuk 1S IS JP 

Apabila koordinat semua titik P(J) dan T(I) diberikan dan semua nilai w(I) juga 
diberikan, maka persamaan (6.108) adalah suatu persamaan simultan yang terdiri dari JP 
persamaan melibatkan JP variabel anu yaitu nilai-nilai yYP(J) untuk 1X J X JP yang tak- 
diketahui. Karena semua koordinat P(J) dan T(I) diketahui maka koefisien pengaruh 
C(I,J) dapat dihitung nilai-nilainya, demikian juga dengan semua nilai-nilai w(T(1),as) 

Ini berarti bahwa nilai-nilai YP(J) dapat dihitung dengan menyelesaikan 
persamaan simultan (6.108) menggunakan salah satu metoda standart penyelesaian 
masalah perkalian matriks, dengan catatan bahwa masing-masing persamaan (6.108) 
(untuk setiap nilai I) itu benar-benar bebas tidak tergantung pada atau dapat ditulis 
sebagai kombinasi linier persamaan lainnya. Dengan perkalian lain nilai-nilai YP(J) yang 
unik dapat dihitung bila matriks koefisien pengaruh C(I,J) tidak bersifat singguler. 

Karena tidak ada persyaratan khusus yang harus dipenuhi dalam penentuan nilai- 


nilai w(I) maka kita boleh saja menentukan bahwa w/(I) bernilai tetap untuk semua nilai 


I. Apabila w(I) bernilai tetap maka semua titik T(I) pasti berada pada suatu garis arus 


tertentu. 

Karena permukaan lingkaran atau aerofoil dan lain sebagainya tidak boleh 
ditembus oleh aliran fluida maka persamaan tersebut sesungguhnya adalah garis arus. Ini 
berarti bahwa bila titik-titik T(I) dipilih berada pada permukaan aerofoil maka 
pembahasan yang telah diberikan sebelumnya dapat diterapkan untuk menentukan nilai- 
nilai YP(J) yang akan memberikan pola aliran sekeliling bentuk aerofoil sembarang. 

Walaupun dalam pembahasan diatas kita tidak menemukan persyaratan tertentu 
yang harus dipenuhi dalam pemilihan panel-panel atau koordinat titik-titik P, Tetapi kita 
tahu bahwa setiap panel singgularitas pasti mengganggu pola aliran angin seragam di 
sekitarnya. Keberadaan aerofoil pasti akan mengganggu pola aliran angin seragam di 
sekitar permukaan aerofoil, tetapi tidak ditempat yang jauh dari aerofoil. Ini berarti 
bahwa panel-panel singgularitas harus berada di dalam atau pada permukaan aerofoil 
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dan tidak boleh berada di luar aerofoil. Suatu pendekatan yang sederhana adalah dengan 
menganggap bahwa panel-panel tersebut berada pada permukaan aerofoil. 

Data tentang geometri atau bentuk aerofoil biasanya diberikan sebagai sejumlah 
terbatas koordinat titik-titik pada permukaan aerofoil tersebut. Titik-titik yang diberikan 
tersebut dapat kita anggap sebagai titik-titik P(J) yang membentuk panel-panel PO(J), 
menghubungkan titik P(J) dan P(J#-1). Panel-panel tersebut membentuk suatu kurva 
kontinyu (mulus) linier sepotong-sepotong yang berbentuk mirip seperti permukaan 
aerofoil yang sebenarnya. Konsep yang diberikan diatas adalah dasar dari metoda panel 
yang akan dibahas lebih rinci dibab berikut. 
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LAMPIRAN 6.1 


Prosedur menggambar kurva iso-nilai fungsi dalam bahasa FORTRAN, dengan semua 
nilai fungsi di titik kisi T(i,j) diketahui. 


Do I-1,IMAX-1 
Do J-1,JMAX-1 
IS —0 


FA - nilai fungsi di titik TA,J) 
FB - nilai fungsi di titik TI,J-1) 
FC - nilai fungsi di titik TI-1,J-1) 
FD - nilai fungsi di titik TI--1,J) 
If (FA"FB.le.0.0) then 
JIS-—1 
XS(IS) XT) 
If (FA.eg.0.0) then 
YSIS)-YTI) 
else if (FB.eg.0.0) then 
YSIS)-YTIH) 
else 
YS(JS) adalah solusi dari G(I,y)-0 di sisi ke 1 
endif 
endif 
If ((FC.EO9.0.0).0r.(FB"FC.It.0.0)) then 
IS-JS-H 
YSS YTI-H) 
if (FC.EO9.0.0) then 
XS(IS)-XTH1) 
else 
XS(JS) adalah solusi dari H(I#1,x)-0 di sisi ke 2 
— endif 
endif 
if (IS.LT.2) then 
If ((FD.E9.0.0).0r.(FD“FA.It.0.0)) then 
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IS-JS-H 
XS(ISYXT(H1) 
if (FD.EO.0.0) then 
YS(IS)-YT) 
else 
YS(JS) adalah solusi dari G(--1,y)-0 di sisi ke 3 
endif 
endif 
if (JS.EO.1) then 
IS-IS-H 
YS(IS)-YT) 
XS(JS) adalah solusi dari H(J,x)-0 di sisi ke 4 
endif 
endif 
if (JS.E0.2) then 
Gambar hasil yaitu garis dari (XS(1), YS(1)) ke (XS(2),YS(2)) 
endif 
enddo 
enddo 
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BAB VII 
METODA PANEL UNTUK ANALISIS KARAKTERISTIK 
AERODINAMIKA BENTUK AEROFOIL 


Dalam bab sebelumnya kita telah membahas bagaimana berbagai macam pola aliran 
dapat diperoleh sebagai hasil kombinasi linear aliran angin seragam dengan panel 
singgularitas tertentu, dimana distribusi kekuatan per satuan panjang singgularitas tersebut 
dianggap diketahui atau diberikan sebagai data masukan. Pembahasan tersebut memang 
bermanfaat untuk memperoleh pengertian dasar mengenai bentuk-bentuk pola aliran yang 
ditimbulkan oleh panel-panel singgularitas seperti sumber, vorteks, dan doblet. Namun 
demikian yang kita inginkan adalah kemampuan menyelesaikan bentuk pola aliran angin 
seragam yang diganggu oleh keberadaan aerofoil. Masalah tersebut dapat dimodelkan 
sebagai angin seragam yang diganggu oleh keberadaan panel-panel singgularitas. Tetapi kita 
tidak tahu pada awalnya bentuk atau distribusi kekuatan singgularitas tersebut. Kita hanya 
tahu bahwa permukaan aerofoil adalah salah satu dari garis-garis arus aliran yang dihasilkan. 
Gangguan keberadaan aerofoil terhadap pola angin seragam dapat diperkirakan hanya 
berpengaruh di kawasan aliran yang tidak terlalu jauh dari permukaan tersebut. Pemikiran ini 
menjadi dasar dari pemilihan panel-panel singgularitas yang harus digunakan dalam 
pemodelan, yaitu bahwa panel-panel tersebut harus membentuk suatu sirkit atau kurva 
tertutup yang tepat sama dengan permukaan aerofoil atau suatu sirkit yang lebih kecil dan 
seluruhnya berada di dalam kawasan yang dibatasi oleh permukaan aerofoil. 

Untuk memudahkan pendekatan penyelesaian masalah, kita dapat memilih panel-panel 
tersebut berbentuk sama dengan permukaan aerofoil yang dipelajari. Data tentang bentuk 
aerofoil biasanya diberikan sebagai sejumlah terbatas koordinat titik-titik pada 
permukaannya. Walaupun jumlah tersebut dapat diperbanyak, misalnya dengan interpolasi 
kubik lentur (lihat bab 4), tetapi bagaimanapun juga jumlah titik permukaan aerofoil yang 
diketahui koordinatnya selalu berjumlah berhingga. Misalnya saja titik-titik PG), untuk 
IKjKITA, adalah titik-titik permukaan aerofoil yang diketahui koordinatnya dan titik P(1) 
berhimpit dengan titik P(JTA). Panel-panel singgularitas yang dipilih adalah yang 
menghubungkan titik-titik Pj), yaitu panel ke j adalah garis yang menghubungkan titik Pj) 
dan P(j#1). Ini berarti bahwa seluruhnya ada JP panel, dimana JP-JTA-1. Selanjutnya kita 
harus memilih jenis singgularitas yang ingin digunakan, misalnya sumber, vorteks atau 


doblet. Sebagai contoh kita dapat menggunakan panel sumber dengan distribusi kekuatan 
seragam pada masing-masing panel, yaitu pada panel ke j kekuatan sumber adalah ()). 
Apabila angin seragam membentuk sudut dengan arah horisontal maka nilai fungsi arus di 
titik T dengan koordinat (XT,YT) adalah 


JP 


W(T) — YT.cosa — XT.sina aliran) 7.1) 
j-1 


dimana C, (T,j) adalah koefisien pengaruh panel kej pada nilai fungsi arus di titik T, yaitu 


OTT, j) 
PT(T,j) 





C,(T,j)--TR(T, Da Jera, i).60(T,j) -PR(T,j).@P(T,j) 
(72) 
poli) - ((xe(j--1)- xe()) #(xe() ve)” 
PO(j).PR(T,j) — (XP(j-t1)- XP(j)).(XT - XP(j)) 
(xe(j--1)- ve(j))(xT — ye(j)) 


PO(i).PR(T, j) — (XP(j--1)— xp(j)).(vT - yP(j)) 
(xe(j--1)- ve(j))(xT -xp(1)) 


OR(T,j) - PR(T,j)-PO(j) 


PTT, j) — (X1 —xp(j)) «v1 20)” 


Va 


OTT, j) — ((xr-— xe(j «7 #(YT- yP(j-1)P) 


0A(T, j) ae Tea) 


OR(T,j) 
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: : 6 Ph PO(j).TR(T, j) 
NT, j) -0@(T,j) -0P(T, j) — tan Pe aa 
(XP(j), YP()) adalah koordinat titik PG). 


Nilai fungsi potensial di titik T adalah 


JP 


HT) — XT.cosa 4 YT.sina-t-— "(i).C,, (Tj) (7.3) 
ja 


dimana GI j) adalah koefisien pengaruh panel kej pada nilai fungsi potensial di titik T, 

yaitu 

C» (Tj) --PO(j)- OR(T, j).InC,,(T,j) #PR(T,j).In PTT, ) #TR(T, j). ACT, j) 
(7.4) 

Nilai komponen kecepatan di titik T, yaitu u yang sejajar arah PO dan v yang searah dengan 

arah normal pada PO diberikan oleh rumus-rumus berikut (lihat bab 5). 





V(T) - va) 4 un u(T,j) — wa 2) 2 


3 1 
A1) Ani) AD 
2x 
Vektor kecepatan Wr, ji) adalah 
Wr,j) —u. tj) 4v.n(j) 
dimana 

li) - cosB(i).xe-sin B(i).y 


EN 


nli) - —sinB(j).xt- cosB(j).y 
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cosB(i) - (XP(j--1)— xP(i))/pO(i) - AxOr(i)/PO(i) 
sinB(j) - (XP(i--1)- vP(1))/po(i) - Ayer(i)/roli) 
Mengingat bahwa vektor kecepatan juga diberikan oleh rumus berikut 
Vj) Xt V,.y 
dimana V, dan V, masing-masing adalah komponen kecepatan searah sumbu-x dan sumbu- 
y, maka 
VE) -u(T, j).cosB(j)  V(T, j).sin B(j) 
NT, j) - VT, j).cosB(j) tW(T, j).sin B(j) 
Sesungguhnya nilai-nilai V, dan V, dapat juga diperoleh dari rumus-rumus berikut 


vr) - Lari) 


- sar (20) - OR(T,j).InOT(T,j) PRT, j).InPT(T,j) Aa ai 


Wlni) Ji) 


- ala) - OR(T,j).InOT(T,j) “PR (T, j).InPT(T,j) erna) 


Kedua cara tersebut akhirnya memberikan hasil berikut 








V,(T,as) K3 coeg) Ara) | 


PT(T,j) 


VI T,as) - ata )).cosp() an Sai | 


Komponen kecepatan angin seragam di titik T adalah 


V,(T,as) - 2 ea(KT.coso 4 YT.sina) 5 cosa 
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V,(T,as) - 2 AK cos #YT.sin a) — sin 
OYT 


Komponen kecepatan gabungan angin seragam dan pengaruh panel-panel singgularitas 


sumber adalah 


(tis ta H-cdr,i) Os) 

V(T) — cosa Sains) (7.6) 
dimana 

ci) ta SA ot) Ati) sing) an 

C,(T,j) - ACT, j).cosB( Dan Ja B3) (9.8) 


Fluida atau udara tidak boleh menembus keluar masuk permukaan aerofoil. Ini berarti 
bahwa permukaan aerofoil adalah sebuah garis arus, jadi nilai fungsi arus di semua titik 
permukaan aerofoil itu sama. Bila titik T dalam persamaan (7.1) ditentukan berada pada 
permukaan aerofoil maka persamaan (7.1) dapat ditulis ulang menjadi 


JP 
»d clj).C, (Dj) -21- we, — 21. (XT.sina — YT. cosa) (7.9) 
-1 


Bila dedak bahwa titik T(i) berhimpit dengan Pi) untuk IXiXJP maka persamaan (7.9) 
dapat diolah menjadi suatu sistem persamaan yang terdiri dari JP persamaan, yaitu 


JP 

rai — 2, —21.(XP/i).sina — YP(i).cosa) (7.10) 

ja 
Persyaratan (7.10) dikenal sebagai syarat batas Dirichlet karena nilai fungsi pada batasan 
medan solusi, yaitu pada permukaan aerofoil, ditentukan bernilai konstanta tertentu 


walaupun nilai w,.,, tidak diketahui. 
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Bila ditentukan bahwa 
Untuk 1«i« JP 


cli,i) - c,(i) a 
xj) “j) | untuk 1xi«s JP 


D(i) — 21(XP/i).sina — YP(i).cosa) 
Cli, JP 41) ——27 
X IP 41) Woo 
maka sistem persamaan (7.10) dapat diringkas menjadi 
JP-4-1 


» Ci i).x(i) Di) untuk 1 six JP (7.11) 


jal 

Sistem persamaan (7.11) melibatkan JP-4#1 variabel anu, yaitu nilai X(j), dan hanya terdiri 
dari JP persamaan, jadi membutuhkan satu persamaan tambahan. Disamping memenuhi 
persyaratan fisis bahwa udara tidak boleh menembus keluar masuk permukaan aerofoil, yang 
dinyatakan sebagai syarat batas Dirichlet (7.11), perumusan masalah juga harus memenuhi 
satu persyaratan aerodinamis yang dikenal sebagai kondisi Kutta. Persyaratan Kutta 
menyatakan bahwa kecepatan aliran di titik tepi buritan atas itu harus sama dengan 
kecepatan aliran di titik tepi buritan bawah, sehingga aliran meninggalkan permukaan 
aerofoil dari tepi buritan secara mulus bersamaan dari permukaan atas dan permukaan 
bawah. Karena panel pertama dan panel terakhir pada umumnya tidak sejajar dan 
membentuk sudut ekor tertentu, maka kondisi Kutta dapat diterjemahkan sebagai 
persyaratan bahwa arah aliran di titik ekor (tepi buritan) adalah searah bisektor sudut ekor. 
Arah garis bisektor sudut ekor dapat dihitung sebagai berikut 


tanB(1)--tanB( JP) 
1-4 tang/(1).tanP( JP) 





tan(B(1) B(3p)) - 
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jadi 





—AXOP/(1).AXOP/(JP) - AYOP/(1).AYOP/( JP) 


dan persamaan untuk garis bisektor adalah 


B(1) 1 B(JP) — an AXOP(1). AYOP(JP) # AYOP/(1). AXOP(JP) | 


y- tang Hp) 190) ).G-x00)) ve) 


Kondisi Kutta kemudian dinyatakan sebagai persyaratan bahwa komponen kecepatan normal 
pada bisektor sudut ekor di titik B, pada bisektor, harus bernilai nol. Titik B dipilih berada 
sedikit di belakang titik ekor, yaitu 


XP — XP/(1)--0.001 
YB -— tan 3 (pt) ar) ot exe) 


dimana (XB, YB) adalah koordinat titik B. 
Bila @B adalah sudut antara garis normal di B dan sumbu-X maka 


OB - 145 (B(1)1BGP)) 


Komponen kecepatan di titik B diberikan oleh rumus berikut (lihat persamaan (7.5) dan 


(7.6) 
V, (B) - cosa ay C,(B,i)-z(i) 


V,(B) - sina - C,(B, j).x(i) 


dan kecepatan normal di titik B adalah 
V,(B) — V,(B).cos0B — V, (B).sinOB 


114 


atau 


(B) —sin(x —0B) lake, i)-cos0B —C,(B,j).sinOB).c(j) 


j- 


V, 


n 


Kondisi Kutta sekarang dapat ditulis sebagai berikut 


Sioarani).xli) — DJP 1) (7.12) 


ja 
dimana 
CIJP-#1,j) — C,(B,j).cos@B —C, (B,j).sinOB untuk 1XSi« JP 
C(JP 41, JP 41) —0 
DJP #1) — —21.sin(x — OB) 
Persamaan (7.11) dan (7.12) kemudian digabung menjadi 
gpu 


»Cij).x(i) —Dli) untuk Isis JP #1 (7.13) 
ja 


Sistem persmaan (7.13) terdiri dari JP#-1 persamaan dan melibatkan JP--1 variabel anu, jadi 
dapat diselesaikan secara simultan untuk menghitung nilai x(j). 

Perlu diperhatikan bahwa gaya angkat suatu bentuk aerofoil itu berbanding lurus dengan 
sirkulasi di sepanjang suatu kurva tertutup yang mengurung bentuk aerofoil. Selanjutnya 
mengingat bahwa singgularitas sumber itu tidak mengimbas adanya sirkulasi, maka model 
yang dijelaskan di atas hanya berlaku untuk kasus aerofoil simetris pada sudut serang « — 0 
yang gaya angkatnya bernilai nol. 

Untuk aerofoil dengan gaya angkat bukannol kita harus menggunakan panel vorteks 


yang mampu mengimbas sirkulasi. Setelah semua nilai c(j) diperoleh maka distribusi 


kecepatan aliran di sepanjang permukaan aerofoil dapat dihitung sebagai berikut, di titik T 
sembarang pada permukaan aerofoil. 


V(T)- Iv (1) stay (9.14) 


dimana V,(T) dan V(T) dihitung dengan rumus (7.5) dan (7.6). Walaupun pada prinsipnya 


titik T dalam rumus (7.14) boleh berada dimanapun juga pada permukaan aerofoil, tetapi 
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